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Abstract 

It is shown that nonlinear terms in equations of gravitons on the back- 
ground of curved space-time of the expanding Universe can solve the prob- 
lem of the negative square of the effective mass formally arising in linear 
approximation for gravitons. Similar to well known spontaneous breaking 
of symmetry in Goldstone model one must take another vacuum so that 
nonzero vacuum expectation value of the quantized graviton field leads to 
change of spectrum for gravitons. There appears two graviton fields, one 
with the positive mass, another with the zero mass. Energy density and 
the density of particles created by gravitation of the expanding Universe 
are calculated for some special cases of the scale factor. Numerical results 
are obtained for the dust universe case. 



Introduction 

There is a problem in quantum theory of gravitons created from vacuum in 
the expanding Universe with nonzero scalar curvature R (inflation, dust, etc.) 
concerning the long wave graviton modes. In linearized theory of quantum 
gravitons in curved space-time of the isotropic homogeneous Universe one ob- 
tains after separation of variables in the wave equation the equation for the 
function dependent only on time. After conformal (Weyl) transformation to 
stationary metric this equation can be understood as equation in flat stationary 
metric with time dependent mass. It occurs that for long waves this effective 
mass squared is negative. All this occurs due to conformal noninvariance of 
the graviton theory for nonzero R leading to tachyonic behaviour of long waves 
modes. The necessity of going from nonstationary metric to the stationary one 
is motivated by finite results for particle creation [T] . In the end one surely must 
return to the original space-time where the obtained results are still finite. 

In some papers (see [5] and references there) it was proposed to consider 
these modes as classical excitations of the field growing in time, so that one must 
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quantize only modes with momentum with the square larger than the negative 
square of the effective mass. However one knows from the quantum field theory 
that tachyonic behaviour disappears if one takes into account nonlinear terms 
neglected in linearized theory. This is typical in theories of spontaneous breaking 
of symmetry due to redefinition of the vacuum leading to its noninvariance to 
this or that transformation of the Lagrangian. In quantum theory based on a 
new vacuum one gets new masses for the redefined quantum field so that there 
is no negative mass square. 

In this paper the analogous program is made for gravitons. It occurs that 
if one is going from the linear theory of gravitons taking into account the next 
order of nonlinearity one gets the redefinition of vacuum solving the problem 
for long wave gravitons. In the result one gets gravitons with zero and positive 
effective mass. 

Differently from the situation in theory of weak interactions where Higgs 
potential with nonlinear term is taken by hand in our case the nonlinear term 
naturally appears as the second order in Einstein equation. 

In the end of the paper the expressions for the particle density and the energy 
density are obtained for gravitons created in expanding Universe with metric 
which has some special dependence of the scale factor on time. 

1 Getting the graviton equation from Einstein 
equation 

Einstein equations in presence of matter have the form 

Rik — ^9ikR = KTik (1) 

or 

K = - \51T) (2) 

Let us consider the case when matter is homogeneous isotropic liquid filling 
the Universe. Then 

Tik = (e + p)uiUk - cjtkP (3) 

where Ui is the four velocity, p - the pressure and s - the energy density of the 
liquid. 

The problem of creation of gravitons in the early Universe was discussed 
in literature with gravitons considered as quantized small term in the metrical 
tensor. Due to absence of exact quantum gravity usually one deals with lin- 
earized analogy with quantization of other quantum fields. First let us obtain 
equations for classical small perturbations of the metrical tensor and then do 
quantization. Consider the graviton pcrturbations-the gravitational waves as 
small term added to the background metric So there are 

(o) 

gik = S tk +hik (4) 
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(°) 

if hik = the 9 ik is the solution of Einstein's equation of the form 

(°) - (°) - 1 - (o) 

R\ = <Tl~-5l T) (5) 
Let us go from the up to low indices and vice verse by using the background 

(o) . (o) . ^ ^ ^ ^ . 

metric 9 ik- h\ — 9 '"/i„fc and expand equations ([2]) in a series in h\: 

(o) . . (°) . 1 . (o) - 1 . 

R I + 6R'k = «( T I - -SI T +5TI - -615T) (6) 
from which due to (O the perturbations hik satisfy equations 

5B}k = - \515T) (7) 

Using the notation (1 + h)^^ for the matrix inverse to (1 + h) with small h], 
(small in the sense that all eigenvalues of the matrix (1 + /i) are smaller than 
the unit)one obtains 

[\ + h)-'{^5l-h}, + KXk~--- (8) 
Write the Ricci tensor and the curvature tensor as 

i?^ = (1 + h)-^\,{-hi ^Ri + \{i + h)-^'M^, + - <f - <f ) 
+1(1 + h)-^\,{i + hrZihikK'''' - iKi - 4n) 

■ iK'''' + 4' - - Knih^' - /.f ^"')))+ R I (9) 

R^k ^R^k +^(1 + h)-'[,{hlk.i + hi.^.i - JV^l^i - h{f,.,) 

■ (/i^; + h^k - hg) - 2hi^Mi - ^u))) (10) 
R =R +(1 + h)'%{-h^^ R7 + {1 + h)~'Ufhf - Kf) 

-/.l,„<^"'+<,(<^"'-/.f^'')) (11) 

(°) 

Here ";" means the covariant derivative in background metric 9 ik- Consi- 
dering in (jlO[) only first degree in hik one obtains the linearized equations for 
hik in (P) as 

(o) (o) (1) 

^In + h;r;k " K^U,^ " /l^;,;„- 5 rk {Jv^ " /^^^i^) + h.k R= ~2k5 T .k (12) 
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Let us consider the background as homogeneous isotropic nonstastionary 
space-time 

ds^ = dt^ - a\t)M (13) 

where 77 is the conformal time. Here the Latin indices take the values 0, 1, 2, 3 
and the Greek - 1, 2, 3. Then write for the scalar curvature and the Ricci tensor 

5R={1 + h)-^l{{l + - ^'\d ^<o,o - %hi, 

<si?^ = (1 + + h)-^\,{Sh};^ + }^^% - <f ) 



+1(1 + h)-'i,ih;i,X'' - (Ki - <n) ■ ihjp' + - ''T 



-2h%^ ■ ih^^ - /if"')) - -^{h^X,^ - 2hl,hi^))) 



2«2^/i.o,o - ^,h\, + ^h$)- ^J^(l + h)--,hio (14) 

where e = ±1, for the closed, open and flat Universe. The sign "f ' is used for 
the covariant derivative in space part of the metric and " , 0" or comma for the 
derivative in conformal time rj. Metric g*'^ is defined up to arbitrary coordinate 
transformations so one can put some auxiliary conditions. 
Solutions ([6|) for hik can be written as 

K^si + v^ + Bi 

where S].,V^,Bl, are irreducible scalar, vector and tensor components of the 
tensor satisfying the conditions 

(°) I (°) (°) i 

g ,k 5'^ = 5^0, g,k = 0, 9 ,k B''' = 

Considering only the gravitational waves exclude the scalar and vector parts 
by putting the gauge conditions 

h = 0, K,^ = (15) 



4 



In this case the Unearized equations for perturbations /i^ take the form 

hlofi + ^hl, + ^h'^,+h'^:; = (16) 

Go from the variables h], to new variables — a{ri)h\ and make the con- 
formal transformation 

9tk ^ gik/a^iv) (17) 

Then one obtains the equation for the field with spin 2 in Minkowsky flat 
space in some external effective field (e = 0) 

/^?,o,o + M^;^-^M^ = (18) 
After separation of variables and Fourier representation ^J.'^{x) 

lilix) = j d^k{g^{r^)e^^^^af, + gl{v)e-^~^'^ af) (19) 

one obtains for the time dependent g^(?7) function the equation 

ffgW + (fc'-^Kr(^)=0 (20) 

which formally has the negative square of the effective mass rri^a^ = —a" /a. 

Differently from the situation in theory of weak interactions where Higgs 
potential with nonlinear term is taken by hand in our case the nonlinear term 
naturally appears as the second order in Einstein equation. 

Equation (|T6|) was considered in the paper of L. P. Grishuk [5]. Calculations 
of gravitational excitations based on eq. (jl8p were made in [5| , where absence 
of the infrared divergence in this method was shown. Calculation of the energy 
density and pressure of created gravitons was made in the papers of A. Starobin- 
sky [5] and V. Sahni [7]. For small k the Fourier transform of the solution (|16p 
was obtained as: 

h-^{k)^a${k) + b0)JJ^^ (21) 

Note however that solutions of the form ([2T|) for small k cannot be inter- 
preted in terms of usual particles, that is why L. P. Grishuk [2] considered them 
as "frozen" modes forming some condensed classical state. 

Here we continue this research considering what changes in the form of the 
condensed state (or the new vacuum) are introduced by next orders in Einstein's 
equations. 

2 Third order equations for gravitational waves 

Let us consider the right hand side of Einstein's equations. For ^ 

5Tl = {5e + 5p)u'uk - 5l5p + {e+p){5u'uk + u'5uk + Su'Suu) (22) 
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(°) 

One has for the four velocities Ui and u i the conditions giku^u = 1 and 
5 ,fc u u " 1. So 

2 fcJu'=+ (5u^(5u''^ (23) 

Note that Sui due to constraints can depend only on hPji'^r^ ... so 5u°'5u^ 
depends on the squares of these terms. But we shall neglect the fourth and 
higher orders. In synhronous reference system vP = l/a,u°' = so from (|23p 

one has 5vP = and 

5Rl ^]^K[5e + iSp), 5R1 = ^K(5^((5p - 5e) 

5R^ ^ an{e + p)6u°' (24) 

Consider the case of flat space e = 0. Then hg~ = hg h^''^ = \h^''^ where 

Greek indices are put up and below by use of the Minkowsky metric. One can 
look on eqs. p4)) as on Euler Lagrange equations for the fields hik- Then due 
to constraints (fT5|) up to terms of the divergence form one can obtain that not 
only /i|rj = but /if /i^rj, = 0, ... so that dH]) can be transformed to 

<5i?^ = (1 + /i)-i"(i(i + h)-^\,(-hgi; - hff) 



-^5^(1 + hr^hi, + (1 + hr^ii^^Ki,,, - ^/i^',o 

-^(1 + M"'!'(/iko<o - HAfi)) (25) 
From (|25p follows that one can take instead of (0 the equations 

{l + h)l5Rf, = \^n{l + h)}{5p~5e) (26) 
Multiply the last equation on 2a^", then from (|25l) . (pS)) one obtains 



+ + h)-'[{hl,hi, - 2/i^o/ij;,o) + (1 + /i)-^;,(/i;:f + h{f) 



^a^K{l + h)'^p{5e-5p) (27) 
Consider first three orders in h\. in equations (|27p . 
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or 



= (5^ + h^^){a^K{5e - 6p) + '^{h\, - hihl)hi^) (28) 

^ a 

+2hi,hr,, + '^Slh-hi^, - '^h^K,) - - 5p) = (29) 

5'pa^K{5p - fe) 

= /^"o/^ko + Kn^t + ^/^r'"/*™,/. + h\.hi; + s'^p^h-x^, (30) 

Putting away the divergence of hl^h\p + \hl'^h\^p + h\,h{f and taking into 

account for fixed nonzero components of the tensor h\. the condition hliK^h^^ = 
after simple transformations one obtains 



So 



",9,0,0 + /3.7 "i'",3,0"i,0 



Now let us go from variables to variables fi^. = a{ri)h\, and make the 
conformal transformation 

Then we obtain the equation in flat Minkowsky space with some effective 
external field 

M/3,0.0 + Mfl'^ f^P S'IA^^O /^r)(M/3.0 

^ ^ ' a a'^ a a 

1 I / ^ l' a.n , l' a n I' CK,n . / ?^,ck\ n /on\ 
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3 Spontaneous breaking of symmetry 
for gravitons 



Let us consider vacuum solution (j32[) depending only on time. In quantum field 
theory this means dependence of vacuum on time. Then 



M?.o,o = -M? - ^^r^U^ (33) 



Taking into account constraints (|15l) in variables /i} = — /i2, M2 ^ /^i o^^*^ 
gets the potential corresponding to (p3)) as 

Write the field close to the minimum of the potential energy 
a 2a4 



- - f4^?(0V^(0), A^^=/i^(0)+?^ (35) 



One must note that the condition (|35p on /^^(O) is the condition of minimal 
energy at some fixed moment to- 

This is the basic idea. Instead of dealing with time dependent m, A we put 
the initial conditions at some to. This corresponds to the principle of minimal 
energy at this moment. Surely m, A at this moment are numbers. 

Take the solution (l35l) as 



/^2(o)-/^?(o) = o, ^,l{o)^-^,l{o)^^,o^\|^ (36) 

Then the Lagrangian 



can be written as 



i = ^(M^(o) + e")'"(/i^(o) + ).n + |^(A^?(o) + ?^)(Mf (0) + ) 

-£(/^?(0)M^(0)+?^Cf + 2/.^(0)Cf)^ 

) - 1^(4/.^ + (c^ef + 4MoC|ef (?! - si) + 8/ig(e! - 



8a 

/I ..2 

'/3 



+4/xgaef +4/ig((c!)^ + (^1)^ - 2e,e2)) 
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Consider the quadratic in terms 



a a 

Taking into account the equality /ig — -r = — and the gauge £} + & — ^ one 

1* a 



obtains the Lagrangian for gravitons 



i(0 = ^C"^L-^((e!r + fe'f) (38) 

which is cahcd by us the effective Lagrangian after the spontaneous breaking of 
symmetry. 

Euler-Lagrange equations have the form 

eJ;„" + ^e^ = o, e2'- + ^e,^ = o, e?- = o, ^^^" = (39) 



One sees that now in (|39| the sign of the mass squared is a correct one. The 
components , are massless while ^ J , ^| have the nonnegative mass squared 
2a" 

. The solutions for diagonal components (or = ^^'^) can be 

a. 

written as ^ the Fourier integral 

/ ^^(Cfc^: {VV^^' + ctg,{n)e-^^^-) (40) 

And one has the equation 

5i: + (fc^ + — )5fc=0 (41) 
a 

This equation is free from the problem of the negative square of the effective 
2a" 

mass if > and one can construct the quantum theory of gravitons based 

a 

on the new vacuum state p6p . One can notice that the vacuum expectation 
value of the field //(?;) is close to the scale factor. For the scale factor a{rj) = rf 
one obtains 

/^o = «(^y)y^ 

One sees that for aU p e [-1;0) (if p G (0; 1) then = -2a" /a > 
and we leave vacuum as fiQ = 0) the dynamical perturbation ^ I which 
contradicts the condition of the expansion of the curvature tensor into a series 
in this perturbation. The value p ~ —1 corresponds to inflation, so we cannot 
deal the inflation model here while other situations can be considered. However 
this case must be considered separately and it is not studied in this paper. 
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4 The Lagrange formalism for gravitons 

One can see from (|40ll4ip that gravitons in the expanding isotropic Universe are 
described by the effective scalar field d{x) with the Lagrangian 

L = V^{'d{x),^^(xy" - ^R'd{x)d{x)) (42) 

where g = det{gik) and R the scalar curvature. There is no factor ^ because 
one deals with two polarizations. Euler-Lagrange equation for the field is = 
^1 = ^ = t9 • a. Look for solutions of (|4T|) in the form pO|) . The numbers c^, cl 
are changed on the operators cr.ci with commutation relations 

k 

[%.ci]=J(fc-^'), [%,?fc,] = [?t,c + ]=0 (43) 
The Fock vacuum state |in > is defined as 

c^rjin >= 0, < in|in >= 1 

Then for d{x) = —£,{x) one has 
a 

k^) = (^27r)Mv) I d''(^k9kiv)e''' + cigk{v)e-'^'n (44) 
where gk{v) satisfy (|4T|) written as 

g'^ + u;lir^)gk^O, ^Hrj) ^ ^ + (45) 



with initial conditions 



1 



9kivo) = , . g'kim) = iV^Mvoj (46) 



The condition for the wronskian 

/ * 



gkMg'k ivo) - glMg'kivo) ^ (47) 

leads to existence of the full set of solutions of (pij) in the sense of the indefinite 
scalar product 

{^i,i2)=i J ^2) (48) 
The Hamiltonian of the quantized field ^(a;) in the metric (|13p has the form 

H{v) = / dx{M' + — (49) 
./ a 



Putting the field ^{x) from into one obtains 
1 

TT^a^iv) Jo 



= ::5747IT / k^dkujk{ri){Ek{ii){ctck +CkSi) 
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+ Fkctc+ + F:ckCk) (50) 
where the coefficients Ek, Ft are expressed through the solutions of the (|45p 

Eu{v) = ^{\9'u\' + \9k?) Fu{7i) = ^{g'^ + gl) (51) 

The corpuscular interpretation can be made in terms of creation and annihi- 
lation operators bk,b^ diagonalizing the Hamiltonian. If 

Ck = al{T])bk- Pkiv)bt 

then the Hamiltonian is 

Hiv) = / k^dkLOkivKEkiv) - mtbk + bkbt) (52) 

« [V) Jo 

The density of created particles and their energy density [T] can be found 
using formulas 

1 

= ^TTT / k'dk\pk\' (53) 
7r^a''(?y) Jo 

5 Some models of graviton creation 

Let us consider some matter filling the Universe with the equation of state p = je 
where p is pressure and e the energy density. One has for the homogeneous 
quasieuclidean isotropic Universe the equation [4] 

— p£ = — then a(r/) = (54) 
Let us take a{r]) = CrfP (take p > 1) and put it into (|45|) . then one obtains 
g"(r;) + (2p(p - 1)-^ + k^)g{in) = (55) 



where 



So the results obtained for the scalar field in [T] are valid for gravitons for 
any scale factor with the scale factor of a given form. In [1] it was shown that for 
the density of created particles and the energy density defined by ([32] -ES]) one 
gets convergent integrals. Let us calculate them. Putting the notation x = krj 
one gets 

^ + (l + ^)ff(^)=0 (56) 
ax X 

Then the energy density of created particles due to ([55]) is calculated as 

£(r;) =< OlTf'lO > 
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1 (l^p + ^\:,)\g[^)\^) - 1) (57) 



The solutions of (1561) arc Bcssel functions 



g{x) = J^J(i - 4m2,a-) + C2xl^Y{\^/l - 4m^,x) 
Then 



« ^TTT^O.OW ^ -— , < m2 < 4 (58) 



1.5 • 10-3i?3/2 

For the density of created particles in the unit volume one gets 



This integral is convergent [T]. For small m (0 < to < 0.5) n(?/) ~ i? and for 
large to (to > 0.5) 71(77) ~ -v/R- 

Consider dust Universe with a{j]) ~ Crf. Then 

4 • 10"^ 

^gr. = ^4 (^0) 

For the background classical matter one has 

_ 2 ■ 10«4 

Smatt. — -j^ (01) 

So for Planckean time {tpi = 10~'*^se/c) the graviton energy density created 
from vacuum is some ten percent of the matter density while at the inflation 
time tinf = 10~^^seA:) it is only 10~^^ of matter. These numbers are consistent 
with our approximation for the metric in perturbation theory. At the modern 
epoch one gets from ([60|) that the energy flow from the time of the end of 

inf 

£ = 0.5-10-^2 ( ) (62) 



inflation i,„f = 10 ^^sek is 



This can be compared with the flow from the Crab nebula [9] . One sees that 
it is much smaller. 

ecrab = lO-' (^^) (63) 
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